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Abstract 

Gibson scaling and related properties of flame-surface geometry in turbulent premixed combustion 
are demonstrated using a novel computational model, Deterministic Turbulent Mixing (DTM). In DTM, 
turbulent advection is represented by a sequence of maps applied to the computational domain. The 
structure of the mapping sequence incorporates pertinent scaling properties of the turbulent cascade. 
Here, combustion in Kolmogorov turbulence (kinetic-energy cascade) and in Bolgiano-Obukhov convec- 
tive turbulence (potential-energy cascade) is simulated. Implications with regard to chemical flames and 
astrophysical (thermonuclear) flames are noted. 



1 



1 INTRODUCTION 

A widely recognized paradigm of turbulent premixed combustion is the propagation of a dynamically 
passive, advected surface at fixed speed mo relative to the fluid. Despite the omission of important flame- 
turbulence interactions such as flame-generated turbulence and stretch-induced variation of the laminar 
flame speed uq, scaling properties deduced from this formulation have proven useful for the interpretation 
of turbulent combustion processes. 

Based on the passive-surface picture, three interrelated scaling properties that describe the structure 
and evolution of turbulent premixed flames (here considering only the flamelet combustion regime) have 
been proposed, and the predictions have been compared to measurements. First, a scale-invariant law 
governing the dependence of ut/uq on m'/mq, where ut is the turbulent burning velocity and u' is 
the root-mean-square turbulent velocity fluctuation, has been derived and confirmed experimentally 
(Pocheau and Queiros-Conde, 1996). In particular, this law conforms to the dimensionally mandated 
scaling ut ~ u' for it' ^ itQ. Second, measurements of the fiame-surface fractal dimension (North and 
Santavicca, 1990) are consistent with the scaling prediction (Kerstein, 1988) D = 7/3. Third, it has been 
proposed that flame-surface fluctuations are suppressed below a length scale Zg at which the generation 
of fluctuations by turbulence and smoothing of fluctuations by flame propagation are in balance. The 
balance scale Ig is determined by setting the eddy velocity equal to uq. The inertial-range scaling 



of the turnover time t{1) of size-Z eddies is invoked, where L is the turbulence integral scale and u is the 
rms velocity fiuctuation. Taking 1/t{1), to be the size-l eddy velocity, /g is determined from l^/rilg) = mq, 
yielding the Gibson scaling (Peters, 1988) 

IJL ~ {ua/u'f. (2) 

Measured flame structure does not obey this scaling. Rather, the available data suggests that lower-cutoff 
behavior is sensitive to flame-turbulence interactions omitted from the passive-surface picture (North and 
Santavicca, 1990; Gulder and Smallwood, 1995; Erard et al, 1996). 

These three scalings are interrelated in that any two of them imply the third, as demonstrated in 



Sec. 3.1. In this regard, the failure of Gibson scaling despite the confirmation of the other scalings has not 
been explained. The approach taken here to address this conundrum is to construct an idealized numerical 
model of dynamically passive, propagating interfaces advected by turbulence, formulated so that all three 
scalings are necessarily obeyed in the limit u' 3> uo- The model is used to investigate the transition to 
the asymptotic behaviors as u'/uo increases. The numerical results suggest an interpretation of the 
experimental results in terms of the influence of flame-turbulence interactions and finite-u'/iio effects. 

The discussion thus far is specific to the kinetic-energy cascade in incompressible, constant-density 
turbulence, governed by the Kolmogorov scaling, Eq. (0). This is the cascade generated by large-scale 
shear applied, e.g., at fiow boundaries. In buoyant turbulent flows with unstable stratification, the large- 
scale energy input is potential energy that is converted to kinetic energy. Potential energy as well as 
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kinetic energy can cascade to smaller scales, and exchange between potential and kinetic energy can 
occur at each scale. It has been proposed (L'vov and Falkovich, 1992) that the potential-energy cascade 
may dominate. This implies the scaling 



consistent with an earlier proposal by Bolgiano and Obukhov in a different context. This 'Bolgiano- 
Obukhov' (BO) scaling is controversial; Borue and Orszag (1997) discuss experimental, computational, 
and analytical results supporting and contradicting Eq. 

Recognizing the considerable if not conclusive support for Eq. (^, numerical results are presented 
for flame propagation in cascades governed by Eq. (^. The inner-cutoff scale and the fractal dimension 
of flames in BO turbulence are predicted from scaling and compared to numerical results. Though this 
regime is of some relevance to chemical combustion, the principal motivation is thermonuclear combustion 
in supernovae. Implications in this regard are noted. 



2.1 General Technique 

A new computational model. Deterministic Turbulent Mixing (DTM), is formulated for the study of mul- 
tiscale phenomena of interfaces propagating in stirred media. Its spirit is to achieve maximum spatial 
resolution by a radical simplification of specific physical mixing processes, e.g., fully developed hydro- 
dynamic turbulence. In particular, mixing is characterized here only by a universal mapping procedure 
for every length scale and by the relative frequency of mapping events. Since interface propagation can 
be implemented as a discrete binary process, a method that fulfills the above description is conveniently 
based on a cellular grid of discrete binary states, henceforth termed 'unburned' and 'burned' in analogy 
with turbulent flame propagation. Furthermore, we have the freedom to choose the spectral distribution 
of mixing events directly and can therefore restrict the calculation to two spatial dimensions without be- 
ing constrained by inverse cascade effects of two-dimensional turbulence. The method is also applicable 
in three dimensions, but this is not implemented here. 

The computational domain is a two-dimensional (A'^ x N) square grid of cells of size (Ax)^, where 
N — 2^ . Every computation starts with the lower half of the domain 'burned' and the upper half 
'unburned'. The algorithm of DTM can then be sketched as follows: during each burning time step 
Tburn , every neighboring fuel cell of a previously burned cell is burned, mimicking laminar interface 
propagation with a speed uq — Aa;/rburn . Here, horizontal, vertical, and diagonal neighbors are included. 
For convenience, we assign Ax = Tbum ~ 1 and therefore uq = 1. Note that only the relative frequency 
of mixing and burning events is dynamically relevant, i.e., the normalization of time and length scales is 
arbitrary. 

Mixing is implemented as a series of discrete events which map the entire grid onto itself. Each 
map Mi convolutes the interface on a characteristic length scale h = 2\ where imin < i < imax < r. A 
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sequence of maps involving k values spanning a large dynamic range leads to a highly complex, multiscale 
interface structure. A time scale Ti = T{li) regulates the frequency of mapping events on the scale h 
relative to the interface propagation time step Tbum • This function relating the time and length scales 
of the problem is the only link to hydrodynamic turbulence; for instance, by choosing a power law 
dependence of n on k we can model different scaling regimes for the turbulent velocities. The goal of 
DTM is to study interface properties subject to various rules for n and rather generic forms for Mi . In 
this approach, complexity arises from a purely deterministic multiscale mixing process coupled to local 
interface propagation. No random elements are involved. 

Owing to the absence of floating-point operations, the computational expense for all the described 
steps is relatively low. Limited mainly by memory considerations, a large dynamic range can be achieved 
and the evolution can easily be followed over many large scale dynamical times ri,„,,^. Identification of 
the largest mixing scale Zmax = 2'"""' with the integral scale L and the minimum scale Zmin = 2'"'" with 
the Kolmogorov dissipation scale Ik yields an estimate for the efi'ective turbulent Reynolds number of the 
Kolmogorov-cascade computations. 



with Ai = imax — imin. Calculations reported here involve Ai = 9, corresponding to Re « 4100 for the 

Kolmogorov cascade. 

On the other hand, this brief description also shows the main restrictions of DTM. First of all, local 
interface propagation is completely independent of the mixing process. Strain effects on mo are omitted. 
Furthermore, feedback of 'burning' onto the Mi and r; is not represented. In particular, there is no 
flamc-gcncratcd turbulence. Most importantly, the maps Mii are purely artificial and are not claimed to 
reproduce detailed mixing dynamics of hydrodynamic turbulence. They can be expected, however, to 
belong to the same universality class as a large ensemble of scale invariant turbulent eddies, and thus 
provide useful insight into the statistical behavior of stirred propagating interfaces. 

2.2 Implementation of Mixing 

Given the general philosophy of DTM as described above, the following mixing properties can be varied 

according to the physical system under consideration: the form of the maps Aii = A4{li), the mapping 
frequencies = T~^(li) (equivalent to choosing a power spectrum of a stochastic mixing process), and 
their normalization with respect to the local interface propagation speed wo- 

We shall flrst describe the procedure for mixing on a single length scale U for a given frequency t^^ . 
The choice of Mi is the least constrained from first principles. Features that should be incorporated into 
discrete maps in reduced turbulence models have been summarized (Kcrstoin, 1991). By construction, 
their effect is intended to be a combined stretching and folding of the interface on the length scale U. 
Consequently, they contain the single parameter U, and their dependence on the spatial coordinates is 
preferably chosen to be continuous and differentiable in order to reflect smooth physical mixing. 
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As a specific representation for Mi , we clioose a sinusoidai deformation of the entire grid, analogous 
to the action of a single Fourier mode. For a given length scale k, the sequence of maps alternates 
between deformation in the x and y directions. The time between successive maps for given li is r^. 
Then for given k, the mapping sequence is: 



for n = 0, 1, . . ., where t denotes the epochs at which the indicated maps are implemented and int[.] stands 
for the nearest integer. Periodic boundary conditions are imposed horizontally, and we assume completely 
burned or unburned material at the bottom and top boundaries, respectively. The computational domain 
is shifted vertically as needed to prevent burned material from contacting the top boundary. 

Figure 1 demonstrates the effect of the alternating x and y-maps, Eq. (^, for a single length scale 
li — r — 3 over a time period of 6 r^. Here, burning is not implemented. Starting with a planar interface, 
the Mi deform it step by step into rolls resembling vortices of a Kelvin-Helmholtz (KH) unstable shear 
layer. Maps acting on multiple length scales can thus be envisioned to produce a superposition of KH 
vortices on all scales. 

This picture is specific to the mathematical definition, Eq. (^), of the mapping sequence. In particular, 
Eq. is formulated so that the point (0, 0) and its periodic images never move. If successive mappings 
are allowed to have different sets of fixed points, more complicated patterns can be obtained for fixed 
h. The intent here is to obtain the simplest possible behavior for fixed k so that multiscale structure 
mainly reflects the assumed distribution of h. For purposes such as the study of flame surface response 
to local turbulence structure, some variant of this formulation might be preferable. A three-dimensional 
formulation can readily be configured to generate helical fluid paths, reminiscent of flow near vortex 
tubes. Any such analogy between the model and properties of turbulent flow would require detailed 
validation that has not yet been undertaken. The intent here is to formulate a model that emulates 
particular scaling properties of propagating surfaces advected by turbulence, irrespective of possible 
broader implications of the model. 

Generalization to multiple length scales is achieved by concurrent implementation of the map se- 
quences A4i{n) consecutively for all i £ [imin,imax]. If events corresponding to more than one scale i 
occur at a given epoch, then the events are implemented in order of increasing i. (Smaller-scale events are 
governed by shorter time scales and therefore are deemed to be 'completed' before concurrent larger-scale 
events. Scaling properties are insensitive to this somewhat arbitrary interpretation.) The dependence 
of the event frequency on k is derived from the notion of a velocity spectrum Vi = v{li) of turbulent 
fluctuations in scale space. We deflne 

(6) 

in analogy with the turnover time of size-Zi eddies. By postulating a power law for the velocity distribution 
in scale space. 




{x,y + mt[li sin{TT x/2li)]} if t — 2nTi 

{x + mt[li sm{Tvy/2li)],y} if t = {2n + l)Ti 



(5) 




(7) 



5 



where A and V are reference length and velocity scales (analogous to the turbulence integral scale and 
rms velocity fluctuation), we can link the inverse mapping frequencies (the time between two consecutive 
mapping events) directly to the length scales: 



Unless otherwise noted, we will henceforth set rj — 1/3, relating the DTM mixing process specifically to 
inertial-range Kolmogorov turbulence. 

In order to normalize the mixing frequencies with respect to the interface propagation, the parameter 



is introduced. As in Eq. (^, the smallest mixing scale is interpreted as the dissipation scale Zk- In 
the derivation of Eq. (^) it was assumed that inertial-range scaling applies at the scale Zg, and thus that 
Ig > Zk- This is equivalent to the requirement z > 1. In all of our calculations, 2 > 2. 

z parameterizes the small-scale mixing. The large-scale mixing is parameterized by u /uq, which is 
defined in the DTM context as u'/wo = L/[uot{L)], where the largest mixing scale h^^^ is interpreted 
as the turbulence integral scale L. This definition omits the contribution of mixing at scales below L. 
A more literal analog of the turbulent velocity fiuctuation u' can be defined within DTM. The choice of 
definition rescales the numerical results, but does not affect the properties of interest here. 

Based on Eq. (^, these parameters obey the relation zu' /uq = ije^^* for the Kolmogorov cascade. 
The turbulent combustion process can thus be parameterized by it'/uo and either Re (the usual choice) 
or z. z is convenient for identifying the range of validity of Gibson scaling, and therefore is adopted here. 
For all Kolmogorov-cascade computations reported here. Re — 2^^, so u /uq = 8/z. 

Common parameters of all computations (except an illustrative case shown in Figure 2) were the grid 
size, A'^ — 2^^ , and the indices imin ~ 1 and imax = 10 determining the scales W = 2'™'" and L = 2*'"'"' of 
the smallest and largest maps, z ranged from 2 to 4, yielding 2 < u' /uo < 4 for the Kolmogorov cascade. 
For the Bolgiano-Obukhov cascade, this z range corresponds to 10.5 < u' /uq < 21. 

Figure 2 shows a typical fiame image from a computation based on the Kolmogorov scaling. The 
multiscale nature of the flame interface is apparent. 



3.1 Kolmogorov Inertial-Range Turbulence 

Here, DTM is used to investigate the interface properties that can be quantified in terms of the box- 
counting statistics described by Menon and Kerstein (1992). This is implemented by partitioning the 
fiame image into square boxes with side lengths Ij =2-',j — 0...r — 2 and counting the number b{lj) of 
boxes that contain at least one newly burned cell. Obviously, b{l()) is the total number of newly burned 
cells, so ut/uq is the ensemble average of b(lo)/N, where b{lo)/N may be interpreted geometrically as the 
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3 RESULTS AND INTERPRETATIONS 
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interface arc length per unit transverse length. Accordingly, Aj = ljb{lj)/{loN) may be interpreted as 
the arc length per unit transverse length of the smoothed interface, for smoothing scale Ij. Henceforth, 
Aj is denoted A{1), omitting the subscript on 

If I is much smaller than the finest scale of surface wrinkling, then the smoothing has no effect, so 
A{1) should converge to a constant value in the limit of small I. The flattening of computed profiles of 
this quantity at small I is seen in Figure 3. 

For I much larger than the flame-brush thickness, the smoothed interface is planar, so A{1) converges 
to unity. The computational domain is not large enough to demonstrate this asymptote, but this is 
immaterial because the unity asymptote is a mathematical necessity. More significant is the tendency 
for the profiles to collapse with respect to l/L with increasing /. This indicates that the flame- brush 
thickness scales as L. 

Interpreting the two-dimensional computation as a planar cut through three-dimensional space, fractal 
scaling implies A{1) ~ where D is the fractal dimension of the flame surface in three-dimensional 

space (North and Santavicca, 1990). This scaling is identifled by linear dependence on a log- log plot of 
A{1). 

The profiles in Figure 3 suggest linear scaling between the inner and outer cutoffs. The linear range 
appears to broaden with increasing m'/mo (corresponding to decreasing z), as one would expect. The 
slope increases with u' /uo, and in fact, the implied fractal dimensions fall within the error bars on a plot 
of measured as a function of u' /uo (North and Santavicca, 1990). The predicted value D = 7/3 is 
not attained by the computed profiles, consistent with the experimental finding that this occurs only for 
u' /uo much higher than the computed cases (North and Santavicca, 1990). 

Next, it is shown that any two of the three scalings discussed in the Introduction imply the third. 
Consider a box-counting profile, plotted in the format of Figure 3, for u' /ito ^ 1. For / <C Zg, it 
approaches a plateau level of order m'/mo because ut ~ u' . For / ^ L, it approaches a plateau level of 
unity. Implicit here is the assumption that the fiame-brush thickness is of order L. However, this is not 
an additional assumption. Rather, this assumption and ut ~ u' both follow from the assumption that 
flame-brush length and time scales are governed by large-eddy length and time scales for u' 3> wo- 

Thus, the ratio of the small-i and large-i levels is w'/wo. In Figure 3, this vertical span is traversed 
over a horizontal span L/lg. If there is fractal scaling, then the proflle has a region of linear dependence 
between these limits. In the limit u' ^ wo, the slope 2 — D of this linear region, based on these estimates, 
is — log('u'/'uo)/ log(L//g). This gives Ig/L ~ {uo/u')^^^^^'^\ showing that the fractal scaling implies the 
Gibson scaling and vice versa. Analogous reasoning shows that if both scalings are obeyed, then ut ^ u' . 

The apparent implication of the computed results in Figure 3 and of the measurements is that fractal 
scaling is more robust than Gibson scaling, because fractal scaling is obtained, but for D values inconsis- 
tent with Gibson scaling. However, a different picture is obtained when Auq/u' is plotted with respect 
a Gibson-normalized horizontal axis, as in Figure 4. Here, /g is defined as {uq/u'Y'L. In this format, 
Gibson scaling, in conjunction with ut ~ w', implies collapse of profiles below the outer cutoff, ut ~ u' 
by itself implies only a collapse of the low-i plateau levels. The profiles indicate some deviation from the 
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latter collapse, suggesting either deviation from this scaling in the u /uq range of the computations, or 
numerical resolution effects. Resolution effects can be discounted only if the computations resolve the 
\ow-l plateau over a significant I range. It is evident that the resolution is marginal by this criterion. 

Whatever the source of low-Z deviation, the collapse between the inner and outer cutoffs indicates 
accurate conformance to Gibson scaling. It is also evident that the universal scaling function thereby 
obtained has significant curvature throughout the range in which it is well resolved. Although an eventual 
transition to fractal scaling with D = 7/3 is a necessary consequence of the convergence of the profiles 
to a universal form (and the assumption that the outer cutoff scales as L) , it appears that the transition 
from the low-Z plateau to the fractal scaling extends over more than three decades of I. 

Thus, Figure 4 indicates that both the Gibson scaling and the relation ut u' are more robust than 
fractal scaling because they are both seen at u /uq for which the latter is not seen. How can this be 
reconciled with Figure 3, and with measurements? 

The apparent linear regions in Figure 3 may be manifestations of the inflection points of the profiles, 
which necessarily occur owing to profile flattening at small and large I. Measured profiles used to infer 
fractal dimensions typically have no greater range of apparent linearity, oftentimes less. Comparison 
of Figures 3 and 4 shows that extraction of fractal dimensions from such data is problematic. A more 
reliable method for determining scaling properties is to demonstrate that a particular scaling collapses a 
family of curves parameterized by u' /uq. 

Onset of fractal scaling may be slow because the mathematical requirement of spatial homogeneity is 
not strictly obeyed. Flame geometry varies with streamwise location. For example, isolated pockets of 
unburned material are more likely to be found on the upstream side of the flame brush (e.g., Figure 2). 
Therefore, flame geometry is approximately fractal only when the flame interface near the flame-brush 
center becomes so wrinkled that it dominates the box-occupancy statistics over a significant dynamic 
range. In this case, local homogeneity near the flame-brush center is sufficient for approximate fractal 
behavior. This condition can be fulfilled only for high u /uq. 

In contrast, Gibson scaling reflects the local balance between wrinkling by advection and smoothing 
by propagation. Because DTM has been formulated so that propagation and advection are both spatially 
homogeneous, this scaling should not vary with streamwise location. This expectation is supported by 
the apparent robustness of Gibson scaling. However, in combustion processes with flame-turbulence in- 
teractions such as flame-generated turbulence, turbulence intensity will, in general, vary with streamwise 
location relative to the flame brush, so the onset of Gibson scaling may be inhibited. 



3.2 Convective Turbulence 

DTM allows freedom of choice for the assumed turbulent velocity spectrum. The impact of a particular 
cascade mechanism on turbulent flame propagation is investigated by fixing the parameter rj in equation 
(0) to the relevant scaling exponent. We used t) = 1/3 in the calculations summarized in Sec. 



3.1 



A particular example of an eddy cascade whose velocity scaling exponent differs from ry = 1/3 is 
buoyancy-driven turbulence. As indicated in Sec. |l| the coupled density and velocity cascades are 
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presumed to be governed by the Bolgiano-Obukhov (BO) scaling 77 = 3/5. The scahng properties of 
premixed flames propagating in buoyancy-driven flows are derived on this basis. Following the arguments 
of Kerstein (1988), we expect a fractal dimension of Dbo = 2|. The lower cutoff scale, denoted iw, is 
the balance scale at which the eddy velocity [l/Lyu' is of order uo, giving 

= {uo/u f'''L . (10) 

DTM provides a simple illustration of the above relations. The results of computations with = 3/5 
are shown in Figure 5. They involved 2; = 2, 3, and 4, corresponding to u'/uo = 21, 14, and 10.5. The 
format of Figure 4 is used, except that / is normalized by Zw in order to demonstrate the BO analog of 
Gibson scaling. Collapse of the computed profiles verifies this scaling. As in Figure 4, it is seen that 
there is a gradual transition to the anticipated fractal scaling. 

These results have signiflcant implications in the context of astrophysics. A specific class of extremely 
bright astronomical events, so-called Type la supernovae, can be consistently explained by thermonuclear 
explosions of compact white dwarf stars (Niemeyer and Woosley, 1997, and references therein). In this 
scenario, the thermonuclear reaction wave propagates from the stellar core outwards as a microscopically 
thin fiame sheet, 'burning' carbon and oxygen to iron-group elements. The resulting stratification of hot, 
burned material underneath the dense outer layers is Raylcigh- Taylor unstable, rapidly giving rise to 
the production of strong turbulence. Hence, the situation fulfills all of the criteria for fiame propagation 
in buoyancy-driven turbulence and, consequently, our results obtained for the BO cascade apply. This 
question is subject to ongoing research; a more detailed analysis will be presented elsewhere (Niemeyer 
and Kerstein, 1997). 

4 CONCLUSIONS 

Previous work has shown that analysis of a dynamically passive, propagating surface advected by tur- 
bulence is a useful framework for interpreting turbulent combustion phenomena. Here, the assumptions 
generally adopted in this framework have been embodied in a computational model. Deterministic Tur- 
bulent Mixing (DTM). 

The model provides an explicit numerical demonstration of the scalings governing turbulent premixed 
fiame structure and burning velocity in the flamelet regime. Moreover this demonstration has several 
implications with regard to turbulent combustion measurements. The main results and conclusions are 
as follows. 

First, the Gibson scaling of the inner cutoff of fiame structure in Kolmogorov turbulence is demon- 
strated by obtaining the collapse of a family of box-occupancy profiles, parameterized by m'/mo, with 
appropriate normalization. This is the first explicit demonstration that Gibson scaling of a propagating 

surface follows directly from Kolmogorov scaling of turbulent advection. 

Second, the structure of the universal scaling function revealed by this collapse is compared to the 
structure of individual profiles. Though individual profiles suggest fractal scaling with fractal dimension 
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D that increases with u' juo, the universal function shows no fractal scaling over the scale range probed 
by the computation. In effect, the different D values of the individual curves reflect the fact that each 
individual curve corresponds to a slightly different range, in the Gibson-normalized coordinate, of the 
universal curve. Thus, the apparent w'/wo dependence of is a measure of the curvature of the universal 
function. 

Third, the results have bearing on the interpretation of turbulent combustion measurements despite 
the simplifications inherent in the model. The results indicate that the measured u'/tto dependence of 
D may in fact be an indication of lack of fractal scaling, rather than variation of a scaling exponent. 
Measurements suggest eventual convergence to D = 7/3, but convincing convergence is not attained over 
a significant dynamic range of u' /uq. Thus, one possible explanation of the empirical failure of Gibson 
scaling is that u' /uq large enough to observe the scaling has not yet boon achieved in the experiments. 
However, model results suggest that Gibson scaling, in the passive-surface framework, is more robust 
than fractal scaling. This suggests that a more likely explanation of Gibson-scaling failure is flame- 
turbulence interactions omitted from the passive-surface formulation. If this is true, then the convergence 
to D = 7/3, if it occurs, may reflect a combination of flame-turbulence interaction effects rather than 
applicability of the scaling analysis. 

Fourth, the model demonstrates the proposed connection between turbulence structure and flame 
structure. When the Bolgiano-Obukhov scaling presumed to govern convective turbulence is adopted 
in lieu of the Kolmogorov scaling, the inner-cutoff scaling changes accordingly. The analysis implies a 
corresponding change of D, but this cannot be demonstrated within the computationally accessible range 
of scales. It is interesting to note that the passive-surface modeling framework may be more relevant 
to buoyancy-driven thermonuclear flames than to chemical flames because the astrophysical flames have 
relatively low thermal expansion, hence weak flame-turbulence interactions (except for gravitational 
forcing), and because m'/wo in these flames is 'astronomical.' 
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Figure 1: Snapshots of the interface during the first 6 iterations of the alternating mapping sequence, Eq. 
(^). For this demonstration, only a single length scale, li — 32, was mapped and interface propagation was 
turned off. The resolution was = 2^. The sense of rotation can be reversed by changing + to — in either 
the top or the bottom line of Eq. 
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Figure 2: Image of a flame in a simulated Kolmogorov cascade after « 30 large scale mapping times. Black 
is burned. The computation used N = 2^, imin = 1, *max = 6, and z = 1. 
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Figure 3: Normalized box occupancy A as a function of length scale, normalized to the integral scale L. The 
symbols represent, from top to bottom: z = 2, 2.5, 3, 3.5, 4. 
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Figure 4: Rescaled box occupancy Auo/u' as a function of length scale, normalized to the Gibson scale Zg. 
The correspondence between symbols and z values is as in figure 3. The line corresponds to a linear decline 
with a slope of —1/3. 
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Figure 5: Rescaled box occupancy Auq/u' in convective turbulence (eddy-velocity scaling exponent r} = 5/3) 
as a function of length scale, normalized to the cutoff scale ly,. The correspondence between symbols and z 
values is as in figure 3. The line corresponds to a linear decline with a slope of —3/5. 
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